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We found an analytical solution for the time dependent filling numbers of the localized electrons 
in a system of two coupled single-level quantum dots (QDs) connected with continuous spectrum 
states in the presence of Coulomb interaction. This solution takes into account correlation functions 
of all orders for the electrons in the QDs by decoupling high order correlations between localized 
and band electrons. 

We demonstrated that several time scales with the strongly different relaxation rates appear in 
the system for a wide range of the Coulomb interaction value. We found that specific non-monotonic 
behavior of charge relaxation in QDs takes place due to Coulomb correlations. 

We also found that besides the usual charge oscillations with the period determined by the de- 
tuning between the QDs energy levels a new effect of period doubling appears in the presence of 
Coulomb interaction at particular range of the system parameters. 

PACS numbers: 73.63. Kv, 72.15.Lh 



I. INTRODUCTION 

The control and manipulation of localized charge in the 
small size systems is one of the most important points 
in nanoelectronicSfii^ Single semiconductor QDs which 
are referred as "artificial" atoms 3 -^ and coupled QDs - 
"artificial" molecules^ are perspective structures that 
may serve for creation of extremely small devices. Sev- 
eral coupled QDs can be used for electronic devices cre- 
ation dealing with quantum kinetics of individual local- 
ized states f^r— Due to this fact the behavior of coupled 
QDs in different configurations is recently under careful 
experimenta l 10 ' 11 and theoretical invest igationi 12 ' 13 

During the last decade vertically aligned QDs have 
been fabricated and widely studied with the great 
success (for example indium arsenide QDs in gallium 
arsenide).— ~— Such experimental realization allows to 
organize strongly interacting QDs system with only one 
of them coupled to the continuous spectrum states. Con- 
sequently vertically aligned QDs give an opportunity to 
analyze non-stationary effects in various charge and spin 
configurations formation in the small size structures^ 

Lateral QDs seems to be better candidates for control- 
lable electronic coupling between two or several QDs by 
applying individual lateral gates. That's why they are 
intensively studied during the last several years both ex- 
perimentally and theoretically 18 ' 19 

Investigation of relaxation processes, non-equilibrium 
charge distribution and non-stationary effects in the elec- 
tron transport through the system of QDs are vital prob- 
lems which should be solved to integrate QDs in small 
quantum circuits.— ~— Electron transport in such sys- 
tems is strongly governed by the Coulomb interaction 
between localized electrons and of course by the ratio 
between the tunneling transfer amplitudes and the QDs 
coupling. Correct interpretation of quantum effects in 



nanoscale systems gives an opportunity to create high 
speed electronic and logic devices i 27 ' 28 In some of the re- 
cent realizations the Coulomb interaction is weak,— but 
for small size QDs the on-site Coulomb repulsion is in 
general strong^ consequently it is important to take it 
into account. In some cases Coulomb correlations can 
determine time-dependent phenomena—i So the problem 
of time evolution of the charge in the coupled QDs con- 
nected with the continuous spectrum states in the pres- 
ence of Coulomb correlations between the localized elec- 
trons is really vital. 

Time evolution of charge states in the semiconductor 
double quantum well in the presence of Coulomb inter- 
action was experimentally studied in.— The authors ma- 
nipulated the localized charge by the initial pulses and 
observed pulse-induced tunneling electrons oscillations. 
Time dependence of the accumulated charge and the tun- 
neling current through the single QD in the presence of 
Coulomb interaction was theoretically analyzed in.— The 
authors described relaxation processes and revealed three 
time rates for localized charge relaxation in the QD cou- 
pled with the thermostat. Several different time rates 
were also found in the system of two and three interact- 
ing QDs coupled with the reservoir <2i~— 

In this paper we consider charge relaxation in the dou- 
ble QDs due to the coupling with the continuous spec- 
trum states. Tunneling from the first QD to the contin- 
uum is possible only through the second dot. We ob- 
tained the closed system of equations for time evolution 
of the localized electrons filling numbers which exactly 
takes into account all order correlation functions for lo- 
calized electrons. It allows to find an exact analytical 
solution for the time dependent filling numbers of the 
electrons by decoupling the high order correlation func- 
tions between conduction electrons in the reservoir (band 
electrons) and electrons localized in the QDs. In such an 
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FIG. 1. Scheme of the proposed model. The system of inter- 
acting QDs is coupled with the continuous spectrum states 
by means of the tunneling rate 7 = Trisot 2 . 

approximation the electrons distribution in the reservoir 
is not influenced by changing of the electronic states in 
the coupled QDs. For QDs weakly coupled to the reser- 
voir the proposed decoupling scheme is a good approxi- 
mation. We found some peculiarities in filling numbers 
for the electrons dynamics arising due to the Coulomb 
correlation effects. 



second dot is empty, then filling numbers for the elec- 
trons in the second QD remain rather small during the 
time evolution of the charge and Coulomb effects in the 
second QD are not so important as in the first one. 

The creation/annihilation of an electron with spin 
a = ±1 within the dot is denoted by cf a /ci a and n a is 
the corresponding filling number operator. The coupling 
between the dots is described by the tunneling transfer 
amplitude T which is considered to be independent of 
momentum and spin. 

The continuous spectrum states are modeled by the 
Hamiltonian: 



pa 

where c£ a /c pa creates/annihilates an electron with 
spin a and momentum p in the lead. The coupling be- 
tween the second dot and the continuous spectrum states 
is described by the Hamiltonian: 



II. MODEL 

We consider a system of coupled QDs with the sin- 
gle particle levels £1 and £2 connected to an electronic 
reservoir (Fig. [IJ. At the initial time two electrons with 
opposite spins are localized in the first QD on the energy 
level £1 (ni CT (0) = n = 1). The second QD with the en- 
ergy level £2 is connected with the continuous spectrum 
states (e p ). Relaxation of the localized charge is governed 
by the Hamiltonian: 

H = H]j + H tun + Hres- (1) 

The Hamiltonian Hd of interacting QDs 

i=l,2cr 

+ £r(c+C2 ff + c lff c+ ), (2) 

a 

contains the spin-degenerate levels Si (indexes i = 1 
and i = 2 correspond to the first and to the second QD) 
and the on-site Coulomb repulsion for the double occupa- 
tion of the first dot. For simplicity we consider Coulomb 
interaction only in the first QD though it is possible to 
obtain closed system of equations for filling numbers cor- 
relators in a general case taking into account Coulomb 
interaction between all the electrons localized in the dots. 
Our model is suitable for the case when the first QD is 
narrow and the second one is rather widei 35 ' 37 Besides, 
if electrons are initially located in the first QD and the 



Hmn = ^2 t(c+ a c 2(7 + c v<y c\ a ), (4) 

pa 

where t is the tunneling amplitude, which we assume to 
be independent on momentum and spin. By considering 
a constant density of states in the reservoir vq, the tunnel 
rate 7 is defined as 7 = TTv^t 2 . 

As we are interested in the specific features of the non- 
stationary time evolution of the initially localized charge 
in the coupled QDs, we'll consider the situation when 
condition (ej — £f)/7 >> 1 is fulfilled. It means that ini- 
tial energy levels are situated well above the Fermi level 
and stationary occupation numbers in the second QD in 
the absence of coupling between the QDs is of the order 
7/ (£2 — £f) << 1 and can be omitted. Consequently the 
Kondo effect is also negligible in the proposed model. 

Our investigations deal with the low temperature 
regime when Fermi level is well defined and the tem- 
perature is much lower than all typical relaxation rates 
in the system. Consequently the distribution function of 
electrons in the leads (band electrons) is a Fermi step. 

We set h = 1 and therefore the kinetic equations for 
bilinear combinations of Heisenberg operators c^/c^ 



cf a c la = nj(t); c+ a c 2a = n%(t); 
cf a c 2 a = n a l2 {t)- c+ a c la = n^(t), (5) 

which describe time evolution of the filling numbers for 
the electrons can be written as: 
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l dt ni 



-T(rf 21 - hl 2 ), 

TfnJ - nj) - (£ + t/«r ff )"2i ~ * 
-T(nl - n?) + (£ + E^Ka - i 7 n? 2 , 



(6) 

where £ = ei — e 2 is the detuning between the en- 
ergy levels in the QDs. The system of Eqs. (|BJ) contain 
expressions for the pair correlators n^ a n 21 and nJ~ CT n^ 2 , 
which also determine relaxation of the localized charge 
and consequently have to be evaluated. In this system 
we neglect high order correlation functions between local- 
ized and continuous spectrum (band) electrons and fulfill 
averaging over electron states in the reservoir. 

Let us introduce the following designation for the pair 

correlators: K a Z ., =< cfci„ct /CV ' > and consider 

only the paramagnetic case < nf >—< n~ a >. Then the 
following relations take place 



One can easily find that Eqs. (|5J) contain expressions 
for the high-order correlators K^x^ anc ^ K%ii\22 ■ 
Their contribution can be easily written in the matrix 
form T: 



T = 
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(12) 
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Since the evolution starts from the initial state with 
two electrons in the first QD and empty second one, 
the system of Eqs. (jHJ for the pair correlators satis- 
fies the initial conditions: #n"n(0) = 1; #2222 (°) = °; 

K a ~, a , (0) = for the other combinations of indexes i, 

131 j 



j. The high-order correlators K1 21 ^ 22 and K 2 ^[ 22 are 
exactly equal to zero due to the fact that they are the 
solution of the linear homogeneous system of equations 
with zero initial conditions. 

The formal solution of the system for the pair correla- 
tors [see Eq. (jHJ] can be written with the help of the evo- 
lution operator. Time evolution of the matrix elements 
Kij [see Eq. ©] is given by the expression: 



*2u? = < ^2>r >=< a» a? >> 



K 



1211 



'12 "1 



'12 



(7) 



The system of equations for pair correlators can be 
written in the compact matrix form (symbol [ ] means 
commutation and symbol { }- anticommutation) : 



Kij(t) = 5>- 45t ) lm iW0)(e l5+ %, (13) 



where H is defined as: H = H + F 
Let us introduce the evolution operator: 



i^- t K=[K,H'] + {K,f} + T, 
where K is the pair correlators matrix 



(8) 



-iHt 



(14) 



Consequently, the time evolution of the pair correlators 
can be found from the following expressions: 
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matrix H has the following form 



H = 
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(9) 



(10) 



-iHt 
„-iHt 



*2m = {^ lHt )i2K{{)) 22 {e tH ^) 22 = $i 2 W*2 2 (t), 



^ = (e-' lHt ) 22 K(0) 22 (e lH ^) 21 = * 2a (t)*ai(t). 

(15) 

Since -^(0)22 in the matrix [see Eq. ©] is equal to 
i^Uii(O) = 1. The evolution operator $ 22 (t) can be 
obtained from the expression for the operator $ 22 (t) by 
the following substitutions: t — > — t and 7 — > —7. Pair 
correlator is a complex conjugate of if 2 m • 

Finally the evolution operators <&ij(t) are determined 
by the equations: 



and the tunneling coupling matrix F is denoted as: 
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|*32(t) 
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with the initial conditions: 



$22 (*) 
$32 (t) 



(16) 
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A. Equations for the time dependent filling 
numbers 

®i j (0) = 6 ij . (17) 

The time dependent filling numbers ni (i) can be found 
The characteristic equation for the evolution operator from the inhomogeneous part of Eqs. @, which results 
&ij(t) eigenvalues Xi has the form: m . 



(H u - X)(H 22 - X)(H 33 - X){H U -X)-T 2 
x[(H u - X)(H 22 - A) + (H u - X)(H 33 - X) 
+(H 33 - A)(ff 4 4 - A) + (H 22 - X){H M - A)] = 0, 



(18) 



where coefficients -Hn, H 22 , H 33 and H^ are deter- 
mined as: 



{[(*|+ n) 2 +7 2 ][(4+^) 2 -e 2 ] 

-4T 2 ( l |-M 7 )> 1 (i) = (*| + 2* 7 ) 
xTU(G^K^ + G?K££), 

where operators G^ 1 and Gj -1 have the form: 



(23) 



Hu = H 4i = -17, 
H 22 = £ + U, 
H 33 = - 2*7. 



(19) 



Each eigenvalue Xi determines the corresponding 
eigenvector: 



(20) 



In our case it is necessary to obtain expressions for 
the evolution operators $12 (t) and $22 (t) with the initial 
conditions $22(0) = 1 and iy(0) = 0. 

Solution for the system of equations which determines 
the functions &i 2 (t) and & 22 (t) can be written as: 



i=l 
4 



-iXit 



-iXit 



(21) 



1=1 



where, constants Cj can be obtained from the initial 
conditions for the system of equations. 



i 


= 0, 


i 


= 1, 


i 


= 0, 




= 0. 



(22) 



(24) 



Solution of the Eq. (|23|) describes localized charge re- 
laxation and consists of the two parts: the first one is 
the general solution of the homogeneous equation (t) 
(right hand part is equal to zero) and the second one is 
the partial solution of the inhomogeneous equation ni (t) . 



ni(t) = nf(i) +n x (t) 



G(t - t )P(t )dt 



(25) 



where G(t — t )- is the Green function of the Eq. (|2"3")l 
with 5(t — t ) in the right hand part, and P(t )-is the 
right hand part of the Eq. (f2"3"|). which appears due to 
the Coulomb correlations. 

General solution of the homogeneous equation has the 
form^ 



n'l(t)=n 1 [A'e- i ^- E ^ t 

+ 2Re(B'e- l{El - E ^ t )+C'e- tiE2 - E ^ t ], (26) 

where coefficients A' , B' and C are determined as: 

\E* -eil 2 



1^2 -Ell 2 



C 



B = - 



2 — -ca| 1-^2 
(Ba-eiX^-ei) 



£ 2 -£i| 2: 



(27) 



1^2 -Sl| 2 

Eigenfrequencies -Ei can be found from the equation: 



{E-e 1 ){E~e 2 + i 1 )-T 2 = 0, (28) 
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FIG. 2. Different time evolution regimes of the filling numbers ni(t) in the first QD in the presence of Coulomb interaction. 
a )- (C + U)h = ((7/7 = 10, C/7 = -10- black line; [7/7 = 5, f/7 = "5- grey dashed line; U/-y = 3, £/7 = -3- black dashed 
line; (7/7 = 1, £/7 = —1- grey line); b). ~ 1 ((7/7 = 10, f/7 = —7- black line; (7/7 = 5, £/7 = —4- grey dashed line; 

(7/7 = 3, 5/7 = —2.5- black dashed line; U/j = 1, f/7 = —0.75- grey line); c). £/7 = ((7/7 = 10- black line; (7/7 = 5- grey 
dashed line; (7/7 = 3- black dashed line; (7/7=1- grey line). Parameters T/7 = 0.6, 7 = 1 are the same for all the figures. 



and have the form 



Ei,2 = -{e x +e 2 -ii) 
±^(ei-e 2 +*7) 2 +4T2. 



(29) 



Green function G(t — t ) of the Eq. (|23|) can be written 



G(f -t') = J2 aie-W-^QQ - i'), (30) 



where A^-are the roots of the characteristic equation 
arising from Eq. (|23|) : 



at 2 +e-i 2 



A1.2 = -ir/ ± 

+ i v / (4T 2 +^ 2 -7 2 ) 2 +4e 2 7 2 ] 1/2 , 
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FIG. 3. Oscillating terms for the two different time evolution 
regimes of the filling numbers m(t). a), f/7 = 0; b). (£ + 
U)h = 0. Parameters r/7 = 0.6, 7=1, U/-y = 10 are the 
same for all the figures. 



- V /(4T2+^_ 7 2 ) 2 +4C 2 7 2 ] l/2 j (31) 



these roots are connected with the eigenfrequencies Ei 
by the relations 



Al,2 — E\ t 2 — E*^ 2 , 
A3 = Ei — E% , 
A4 = E2 — El . 



(32) 
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(33) 



Coefficients a; are determined as: 



Let us now focus on the two limit cases when the 
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expressions which determine the dynamics of the filling 
numbers have a rather compact form. The first one cor- 
responds to the situation when the detuning between the 
empty energy levels in the QDs is equal to zero: £/7 -C f . 
The second one deals with the situation when the sum of 
the detuning and the half value of Coulomb interaction 
is equal to zero. This means that the resonance between 
the half occupied energy level in the first QD and the 
empty level in the second QD takes place: {£+U)/f <C 1. 
We shall also consider that in both cases the condition 
T < 7 < U is fulfilled. 



For rj <c 1 the time evolution of the filling numbers 
n\(t) can be determined by: 



ni(*) 



1 



. T 2 „ , T ! T 2 

, ^ (e -^-e-^) + 0(^). (39) 



c K + u)h « i 

In the second case of interest ((£ + J7)/7 <C 1 but 
P/7> 1) eigenvalues are: 



B. f/ 7 < 1 

The eigenvalues of the characteristic equation in the 
first case (£/7 -C 1) within the accuracy ^ have the 
form: 



Ai = U-i 

A 2 = -27 - i 



2T 2 7 



C/2 ' 

.2T 2 

1 

7 

2T 2 



A3 = — 2«7 — i , 

7 

A 4 = -17. 

So, the evolution operators can be written as: 



(34) 



$12 (f) = — (e - e ' 



*»(*) = (1 



2T 2 



)> 
2T 2 



)e ^ + -^-e 



(35) 

and time dependence of the pair correlators i^lm an{ l 
i^X2ii is determined by the product: 

«2in(*)=*ia(*)*52(*), 
«Su(<) = W (36) 
Expression for the P(t) [see Eq. ([23)1 ] in the case of the 
resonance between empty levels £,/"/ — has the form: 

P(t) = 4T 2 7e~^* + 2T 2 Ue-^ cos(Ut), (37) 

2 

where 77 = jj^p-s ■ F° r V — 1/2 the inhomogeneous 
part of the time evolution of the filling numbers Hi(i) 
can be written as: 



ni(t) 



T 2 



[(-27* 



-* +e -27* 



4e -^* _ 4e-7*] + ±L_1^ — [cos(Z7t) - 1] 



0( 



U 2 U' 



(38) 



Ai = 


U 

2 "* 


8T 2 7 

f/2 > 




A 2 = 


—17 + 


8T 2 
17 ' 




A 3 = 


-217- 


J7 
f 2" 


i8T 2 7 

+ 

t/ 2 


A 4 = 


—17. 







4T 2 
[7 ' 



Evolution operators have the following form: 



(40) 



*i2(i) = jj{e 2 
$ 22 W = (l-§J)e-^- 



(41) 



When the condition (£ + U)/j = is fulfilled, P(i) 

within the accuracy ^5 and ^ is determined by the ex- 
pression: 



P(t) = -T 2 C7 2 (e^*-^ + /i.e.) - 4T 2 7 V 



(42) 



The inhomogeneous part of the time evolution of the 
filling numbers n\ (t) with the accuracy jp has the form: 



ni(f) = --(l 



-14 



772 \ e jja 



t 2 r/ r 2 

-2 w e^si n( -t) + ( ^). 



(43) 



It is necessary to point out that relaxation of the 
filling numbers in the proposed model can be ana- 
lyzed by means of more simple method — the self- 
consistent mean-field approximation^^ In this approx- 
imation correlation functions U (n~° nfj) in the Eqs. © 
are substituted by the expressions U(n^ cr )(nfj). Such 
substitution is valid in the case when filling numbers for 
the localized electrons n~ a change their values rather 
slow. Calculation scheme consists of the two steps. On 
the first step one has to substitute the initial energy level 
position £j by the expression ei = £j + U (n~ a ) and to 
evaluate the time dependent filling numbers. The second 
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FIG. 4. Time dependence of the filling numbers for the elec- 
trons n\{t) in the presence of Coulomb interaction: compar- 
ison of the exact solution and the mean-field approximation. 
Black line corresponds to the exact solution, black dashed 
line corresponds to the mean-field approximation. Grey line 
demonstrates relaxation of the localized charge in the absence 
of Coulomb interaction, a). U /y = 5, f/7 = —3; b). U/y = 3, 
£/7 = —2. Parameters T/-y = 0.6, 7=1 are the same for all 
the figures. 



step deals with the self-consistent calculation of the time 
dependent filling numbers for the electrons. For some 
ranges of the system parameters mean-field approxima- 
tion reveals qualitatively good results. 35 But in general 
case the mean-field approximation is insufficient to de- 
scribe the relaxation processes in the system with corre- 
lations. 



III. RESULTS AND DISCUSSION 

Time evolution of the filling numbers for the electrons 
strongly depends on the relations between the system 
parameters. 

If condition K + ^)/7<l is fulfilled, the Coulomb in- 
teraction value increasing leads to the decreasing of the 
filling numbers relaxation rate [see Fig. HJa)]. For the 
large U relaxation rate is rather slow and is of the order 

of 7nonres = which is typical for the system of two 

coupled QDs without Coulomb interaction with |£| ~ U. 
By the decreasing of the Coulomb interaction value U 
we achieve the situation of resonant tunneling between 
the localized states and consequently relaxation rate be- 
comes larger. On the Fig. |2fc) the situation of resonant 
tunneling between the empty energy levels = is 
demonstrated. In this case the relaxation of the local- 
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FIG. 5. Relaxation of the filling numbers n\{i) in the presence 
of Coulomb interaction in the case of the resonant tunneling 
between the empty energy levels in the QDs. Black line cor- 
responds to the exact solution, black dashed line corresponds 
to the mean-field approximation. Grey line demonstrates re- 
laxation of the localized charge in the absence of Coulomb 
interaction, a). 5/7 = 0, U/7 = 3; b). 5/7 = 0, U/j = 1. 
Parameters T/7 = 0.6, 7 = I are the same for all the figures. 

ized charge takes place with the typical rate very close to 
the value 7 res = 2— and is almost independent on the 
Coulomb interaction value. Let us notice that relaxation 
processes are governed not only by the typical exponents 

e 7 and e 1 but also by the pre-exponential factor, 
which linearly increases in time in the resonant case [see 
Eq. ((Ml)]. 

A very special relaxation regime exists in the system 
if condition ~ 1 takes place [see Fig. HJb)]. In 

this regime Coulomb correlations result in formation of 
a dip in the time evolution of the localized charge. At 
the initial relaxation stage the charge in the first QD 
rapidly decreases due to the almost resonant relation be- 
tween the level in the second QD and effective single elec- 
tron energy in the first dot. It follows from the third 
and the fourth Eqs. ^ of the system that changing 
of the effective energy levels detuning is determined by 

U Rel^yS^M^-] which differs from the typical mean- 

field expression U (h~° '(f))j2£ 

At a certain instant of time the effective single elec- 
tron level falls down beneath the level in the second 
QD. At this moment the inverse charge begins to flow 
from the second QD to the first one. The occupation 
in the first QD demonstrates significant increasing after 
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reaching minima value (the dip formation) . Filling num- 
bers almost reach the initial value for the large values of 
Coulomb interaction. After the dip formation the typical 
time scale which determines relaxation of the filling num- 
bers is close enough to the value 7„ nres = 2-^^. This ex- 
planation gives qualitative picture of the dips formation. 
The exact solution shows, that Coulomb correlations are 
responsible for such non-monotonic behavior. This effect 
is determined by the inhomogeneous part of the exact 
solution for time evolution of the filling numbers in the 
first QD [see the first term in Eq. (|43[)] ■ And this in- 
homogeneous part appears due to complete account for 
time dependence of the high order correlators [P(t) in Eq 
(1231) and Eq. (US]))]. That is why time evolution of the 
filling numbers for the electrons differs considerably from 
mean-field approximation. The width of the dip can be 
roughly estimated as 1/8 • 7^ 1 „ T . es ■ 

We would like to stress that the non-monotonic be- 
havior, which we discussed above, is not connected with 
the usual quantum oscillations between two energy lev- 
els. Such oscillations also take place during time evo- 
lution, but the amplitude of these oscillations is rather 
small (of the order jp ) . Only these small oscillating con- 
tributions to the total electron density are shown on the 
Fig. [3] Oscillations are always present in the case of 
strong Coulomb interaction for all the values of the ratio 
T/j. We found out that besides the oscillations gov- 
erned by the system parameters T and U, oscillations 
with the double period exist in the system. Oscillation 
period doubling is mostly pronounced in the case when 
resonant tunneling takes place between the half occu- 
pied energy level in the first QD with the initial charge 
and empty level in the dot coupled with the continuous 
spectrum states (£ + U)/-f <C 1 [see Fig. E^b)]. Double 
period oscillations disappears with the decreasing of en- 
ergy levels detuning £. In this case oscillations period is 
determined by the value of the Coulomb interaction [see 
Eq. ([Mill- 
Comparison between the exact solution and the mean- 
field approximation is demonstrated on the Figs. 0HSJ It 
is clearly evident that both methods reveal such similar 
peculiarities of the system behavior as several time ranges 
with considerably different relaxation rates. For some 
ranges of the system parameters formation of the dip can 
be also reproduced in the mean-field approximation (see 
Fig. S|) . Figure 0] also demonstrates similar behavior of 
the exact and the mean-field solutions at the initial stage 
of relaxation. But the dip reproduces incorrectly in the 
mean-field approximation. 

In the case of resonant tunneling between the energy 
levels in the QDs (£/7 = 0) the exact solution and the 
mean-field approximation reveal strong mismatch [see 
Fig. 03a)]. Exact solution demonstrates rather smooth 
time evolution of the localized charge while the solution 
obtained by means of the mean-field approximation re- 
veals abrupt changing of the localized charge amplitude. 



If the Coulomb repulsion decreases the correspondence 
between the exact and the mean-field solutions becomes 
better [see Fig. EJb)] 

Finally let us return to the influence of the Coulomb 
repulsion in the second QD on the evolution of the filling 
numbers. In this situation time evolution of the filling 
numbers for the electrons can be analyzed by means of 
the equations obtained for the model when Coulomb in- 
teraction acts in the first QD [see Eq. ©] after substi- 
tuting the value Un^ a by the Uh^ in the Eq. @. The 
results are shown on the Fig. [6] and it is clearly evident 
that in this case the influence of Coulomb correlations on 
the relaxation of the filling numbers is rather weak. 

IV. CONCLUSIONS 

We have studied time evolution of the filling numbers 
in the system of two interacting QDs coupled with the 
continuous spectrum states in the presence of Coulomb 
interaction in one of the dots for a wide range of the 
system parameters. The solution describing the system 
dynamics was analyzed in the assumption that the band 
and localized filling numbers for the electrons are uncou- 
pled. This solution exactly takes into account all order 
correlators for the localized electrons in the QDs. 

We found strongly different relaxation regimes in the 
system of coupled QDs depending on the ratios be- 
tween the system parameters. Interesting manifestation 
of Coulomb correlations is the formation of the dip in the 
time evolution of the localized charge. Such reentrant 
charge behavior is not the result of simple quantum os- 
cillations between the two energy levels. Oscillations of 
this type are also present in the system but have small 
amplitude in the case of the strong Coulomb interaction. 
Interaction effects lead to the appearance of oscillations 
with double period at particular range of parameters to- 
gether with the oscillations governed by the detuning be- 
tween the energy levels. 

We compared our results with the mean-field approxi- 
mation. The mean-field approximation can give in some 
cases qualitatively similar peculiarities of the system be- 
havior: several time ranges with considerably different re- 
laxation rates and dip's formation. But in many regimes 
the results of the mean-field approximation do not coin- 
cide with the exact solution. Even if the mean-field ap- 
proximation qualitatively correctly predicts appearance 
of the dip, it's shape and width strongly differs from the 
exact solution. 
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FIG. 6. The influence of the Coulomb interaction in the second QD on the time evolution of the filling numbers n\(t) in the 
first QD. Black line corresponds to the exact solution. Grey line demonstrates relaxation of the localized charge in the absence 
of Coulomb interaction. a).£/7 = —5; b).£/7 = —4; cj.^/y = 0. Parameters T/y — 0.6, 7=1, U/y = 5 are the same for all the 
figures. 
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